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. We study the control by electromagnetic fields of molecular alignment and orientation, in a linear, 

' rigid rotor model. With the help of a monotonically convergent algorithm, we find that the optimal 

field is in the microwave part of the spectrum and acts by resonantly exciting the rotation of the 
. molecule progressively from the ground state, i.e., by rotational ladder climbing. This mechanism is 

' present not only when maximizing orientation or alignment, but also when using prescribed target 

states that simultaneously optimize the efficiency of orient at ion/ alignment and its duration. The 
extension of the optimization method to consider a finite rotational temperature is also presented. 



O 



I. INTRODUCTION 



C/3 

O ' External fields can be used to manipulate molecules to achieve molecular axis alignment or orientation. Here, 
' ^ alignment refers to setting the molecular axis parallel to a laboratory fixed frame, while orientation implies that 
the molecular axis has in addition the same direction as the laboratory fixed frame. These two goals have a wide 
range of applications in fields such as chemical reactivityji surface processingy^'^ nanoscale design;^ attosecond pulse 
production, ^'^ and quantum information processing^ 

Efficient alignmen1i2ii£ and orientationiiiiS can be achieved by laser-induced adiabatic passage from field-free rota- 
fS| ' tional states to aligned pendular stateSfiii^ji^^iSiii but it is lost at the end of the laser pulse. Field-free alignment 
J> is possible by sudden excitation using pulses much shorter than the rotational period of the molecule. — 
On ' Achieving orientation with short pulses is more difficult since spatial symmetry breaking is required to give the direc- 
\ tion of orientation, but it can be done using half-cycle pulsesSSiSi or specially tailored laser pulsesiSS*2& This has led 
to different proposals for alignment and orientation using series of short impulses (kicks) iSLSSiSiSSiSiiSSiSSiM See also 
the recent review of the subject by Stapelfeldt and Seideman, and references thereini2^ 
{p^ ' The purpose of the present study is to find the electromagnetic fields that produce the best possible orientation 
\ or alignment. We start by presenting in Sec. ^ the rigid rotor model used to describe the rotation of a linear 

■ molecule, along with the cost functionals that describes the required control objectives, in terms of both observables 
\ measuring orientation or alignment and target states that embody the efficiency of orientation/alignment along with 
. its persistence. 

J>-j' The optimization procedure itself is based on monotonically convergent algorithms^Si2£i2& that are guaranteed to 
i-S^ improve at each step the cost functional chosen. The corresponding algorithm for the control of alignment/orientation 
is presented in Sec. Ill CI 

^ , As we will see in Sec. IIIII the fields leading to optimal orientation and alignment are in the microwave part of 
k>( ' spectrum and lead to rotational ladder climbing, i.e., the molecule is resonantly excited successively from one 

r> \ rotational level to the next. The possibility of controlling rotational excitation by ladder climbing using microwave 

■ fields was first proposed by Judson et alm^^ A process similar but resting on Raman excitation of ro-vibrational 
" " states with chirped pulses has been used to create an optical centrifuge for molecules4ii^242iilii^ 



II. MODEL 



A. Time-dependent Schrodinger equation 



The dynamics of the molecule interacting with the electromagnetic field is obtained by solving the time-dependent 
Schrodinger equation (TDSE). We restrict ourselves to the case of a linear molecule in a rigid rotor approximation. 



2 



yielding the Hamiltonian (in atomic units, h = 1) 

H = BJ^ - fioSit) cose - [{an - a±) cos^ 9 + aj] (1) 

where B is the rotational constant, J is the angular momentum operator, 9 is the polar angle positioning the molecular 
axis with respect to the polarization vector of the linearly polarized electric field of amplitude £{t), /io is the permanent 
dipole moment, and a|| and a± are the dipole polarizability components parallel and perpendicular to the molecular 
axis, respectively. Because of the cylindrical symmetry about the field polarization axis, the motion associated with 
the azimuthal angle can be separated and M , the projection of the total angular momentum J on the axis, is a good 
quantum number (AM = 0). The TDSE is solved numerically starting from the ground rotational (isotropic) 
state J = M = 0, using a basis set expansion of the wave function -0 in terms of the spherical harmonics Yj^m, 

oo 

^{0,t)^J2''j(*^^J'O{0), (2) 

J=0 

the cj being complex coefficients and the coupling terms due to /io f^nd a being then analytical. For computational 
purposes, only the first 10 terms in the sum in Eq. jSJ are kept, and we have checked that the results arc not affected 
by using a bigger basis. 

Because of the presence of both the dipole moment /i and the polarizability anisotropy Aa = q;[| — a±, the results 
are not molecule-independent. However, the role of the polarizability is negligible for the fields obtained, making 
them applicable to any linear molecule with a proper scaling of the amplitude and frequency of the electric field. The 
results will thus be presented with time expressed in units of the rotational period Trot = h/2B, the electrical field 
as hqS/B, and energy as E/B. The parameters actually used in the calculations are those for the HCN molecule: 
B = 6.6376 X 10"*^, n = 1.1413, ay = 20.055, and a± = 8.638 (all in atomic units). 



B. Cost functional 



As we are seeking to optimize molecular orientation or alignment, our cost functionals are based on their respective 
measure, the expectation values {cosd) and (cos^0). A molecule will be oriented when |(cos0)| ~ 1, with the sign 
indicating in which direction it is pointing; an angular distribution symmetric with respect to 9 = 7^/2 will yield a 
value of zero. The expectation value of cos^ 6* is 1 when the molecule is aligned, starting from 1/3 for the isotropic 
case. 

The first case we consider is a cost functional of the form 

MS) = mf)\ 6 10(if)) - t \{t)S\^) dt, (3) 

Jo 

with O an operator chosen to be cos 9 + 1 for orientation and cos^ 9 + 1 for alignment, the identity operator / being 
used for convenience (e.g., it ensures that O is positive) without modifying the extrema of J7i, and t[ the time at 
which the interaction with the field ends. The last term in Eq. © is a penalization on the amplitude of the field, with 

This imposes a strong constraint on the maximum amplitude of the electric field, such that the field strengths are 
comparable to those that can be achieved for half-cycle pulses^^i^ allowing a comparison with previously published 
results for alignment and orientation using kicksiS2iS^i2L22422i2242ii22iSi^ In addition, the form of Eq. ensures a 
smoother, and thus more realistic, turn-on and turn-off of the field. We point out that, in any case, this penalty term 
is not an essential ingredient of the monotonic algorithm and can be relaxed (e.g., to allow fluences achievable with 
laser pulses) or completely eliminated. 

The downside of such a cost functional is that it takes into account only the efficiency of the orientation/alignment, 
not its persistence. Once the field is turned off, the free rotation of the molecule will lead to the disappearance of the 
orientation/alignment as the different J components in the wave function dephase, followed by revivals at intervals 
of one rotational periodii^ The best orientation/alignment is obtained by confining the molecule to a narrow angular 
distribution A9, which corresponds to exciting a broad rotational band A J by referring to an uncertainty principle 
A J • A9 ~ The problem is then that, conversely, a broad rotational spectrum exhibits narrow features in the time 
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domain, i.e., the greater the orientation/aUgnment, the shorter its duration. A compromise has thus to be made, as 
can be achieved by considering the best orientation/ahgnment possible for a restricted maximum r otation al excitation. 
The procedure on how states with such characteristics can be obtained is given in detail in Refs. where it can 

also be seen that Jmax = 4 offers a good compromise, leading to an orientation of (cos0) « 0.906 or an alignment of 
(cos^ 9) « 0.837, both lasting of the order of 1/lOth of the rotational period. The cost functional is now 

J2 {S) = 23? (V'target | V'(if )) ' f ^ A(t)f ^ (i) dt, (5) 

JO 

where V'targot denotes the target state corresponding to orientation or alignment, as given in Tabled and K the real 
part. Note that because of the norm conservation properties of the Schrodinger equation, the cost functional © has 
the same minima and critical points as 

J{£) = - II^Atargct - V'(if)f - / * ms'it) dt, (6) 

Jo 

which measures the distance between V'targot and ip(t{). 

In all cases, the time at which the field is turned off and the cost functional measured is chosen as t[ = 9.5 x 10^ a.u. w 
20ri.ot for the results presented here. Shorter durations lead to results either similar or less significant. 



C. Monotonically convergent algorithm 



The algorithm used to find the optimal field is based on a general class of monotonically co nverge nt algorithms 
recently proposedi^ We present here the algorithm associated to Ji and refer the reader to Refs. ISOlSlI for a detailed 
discussion the algorithm in a time-discretized framework. At the maximum of the cost functional J7i, the Euler- 
Lagrange critical point equations are satisfied; a standard way to write these equations is to use a Lagrange multiplier 
x{0,t) called adjoint state. The following critical point equations are thus obtained: 



^dtx = Hx, xiti) - 0(^(tf)), 

X{t)£{t) = -3? (x(0 1 MO cos 61 + 2£:(t) (Aacos^ 6l + a±) IV-W) , 



(7) 



where 3 is the imaginary part of a complex number and tpo the initial state of the controlled system. 

Given two fields £ and £' and the corresponding states ip, ip' and adjoint states Xj x' defined by Eq. 10, one can 
write 

AJi = J,{£')~M£) 

= (V/(tf)-^(tf)|0|^'(tf)-VXif)) 

[£'{t) - £{t)] I 23 {^'{t)\ ^0 cos \x{t)) 



+ [£'it)+£{t)] 



23 (^'(i) 



Aa cos^ 9 + a± 



\xit)) - x{t) 



dt. 



(8) 



The first term of this sum is positive since both choices O = cos+/ or O — cos^ +/ are positive. Given £, the 
integrand provides thus an implicit criterion in terms of £', the satisfaction of which guarantees the positivity of AjTi. 
An explicit choice of £' can be exhibited: the integrand of Eq. JHJ is a second-order polynomial with respect to £' and 
for a large enough value of A(i) the coefficient 23 ^"^°"^^+°^ \x{t)) - X{t) of £''^{t) is negative. It has thus a 

unique maximum, given by the cancellation of the derivative. The value obtained by this method is 

3(^'(t)|^ocos0|xW) 



£'{t) = - 



23(vywi 



Aa cos^ 9'\-a± 
2 



xit))-m 



(9) 



The algorithm derived from the previous computations is then given by the following procedure: given at step k a 
field £'^ and its associated state tp'^ and adjoint state x*^, compute simultaneously 5*^+^, tp'^'^^ by 



j(^'= + Ht)|M0 cose|x'=(t)) 



idti^''+'^{t) 
^fc+i(6»,i = 0) 



B - fi„£''+^{t) cost 

MO)- 



\x''(t))-m 



(Aacos2 6i + a^) t/j'^+^t), 



(10) 
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Then compute backward evolution of x'^'^^ by 

i idtx''+Ht) = [b - fio£''+Kt)cose- (Aacos^g + a^)] 

The arguments above show that 



(11) 



ji{£'+')>jii£n- (12) 

III. RESULTS 
A. Optimizing orientation 

The electric field obtained using the cost functional Ji [Eq. ^] with the observable O — cos 9 + I, i.e., for the 
optimization of the orientation, is given in Fig. ^a). To better analyze the result, we have performed a short-time 
Fourier transform (STFT),S 

p+oc 

J^{i^,t)= £{t)w{t ~t)e-''^'"'^dT, (13) 



where w is a Tukey-Hanning window with a temporal width of 1.9 x 10^ a.u. The frequency distribution J- can be 
seen in Fig. ^b), where the ordinate is the dimensionless value 2t'Trot, corresponding to the dimensionless energy 
E/B. The energy spacing between rotational states J and J + 1 being 25 J, we see clearly from Fig. ^b) that the 
field is initially resonant with the J = —> 1 transition, and subsequently comes in resonance with higher and higher 
pairs of rotational levels J=1^2, J = 2^3, ... Looking at the population of the rotational states, in Fig. |2 
we indeed find that, starting from the ground state J = 0, the molecule is pumped to the first excited state, then 
to the second, etc. At the end of interaction with the field, the population distribution is such that an orientation 
of {cosO) {t{) — 0.909 is attained (Fig. |3Jl. In other words, the molecule is oriented by a process of rotational ladder 
climbing. 

If we instead take the cost functional J2 [Eq. Q] with V'targot the target state for orientation given in Table 
we obtain a result very similar to the previous one, as shown in Fig. 21 The main difference is the absence of the 
frequency component at 2i/Ti.ot = 10, which is easily understood from the fact that it corresponds to the J — 4 ^ 5 
transition, while the target is restricted to Jmax = 4. The resulting dynamics of (cos 9) are nearly indistinguishable, 
as is seen in Fig. |21 The similarity of both results can also be explained by looking at the projection on the target 
P = I ("(Atargct |V'(^f )) 1^- For the wave function obtained for the optimization of the observable, we already have 
P = 0.9933, the optimization of the projection on the target allowing to reach P — 0.9969. This efficiency is better 
than that obtained when kicking the molecule with short pulsesj^^iS but the time necessary to reach the optimized 
state is of the order of 20 rotational periods, while it takes less than one molecular rotation with kicks. The time 
needed for ladder climbing also explains why the state reached for the optimization of {cos 6) is almost the same as 
the target state, since the time limit imposed by tf constrains the maximum value of J than can be excited. It can 
be likened to the reduced Hilbert space used when defining the target statei^^*^ 



B. Optimizing alignment 

The result for the maximization of the operator O — cos^ 9 + I to achieve alignment is given in Fig. El The field 
obtained is almost the same as the one obtained for orientation, except that the rotational excitation happens at a 
quicker pace, as displayed in Fig. \7\ where it seen that J = 5 is now significantly populated. The alignment obtained 
is (cos^ 9) = 0.866. 

Changing now the the target state for alignment (Table P) , we see in Fig. |S1 that the field obtained is significantly 
different. The frequency component corresponding to an energy of E = AB is present for a longer time and components 
ai E = 6B and, to a lesser extent, E = 2B, reappear near the turn-off of the field. The time dependence of the 
population of the rotational states in Fig. (51 gives the explanation of this phenomenon: the populations of J = 1 and 
3 are pumped down by these later components, since only even J levels are populated in the optimally aligned target 
state. The original excitation to the odd levels was necessary as the rotational states are only significantly coupled 
via the permanent dipolc moment, implying the selection rule AJ = ±1, the role of the polarizability being here 
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negligible. This excitation-deexcitation scheme leads to a projection on the target state of P = 0.9950, compared 
with P ~ 0.5487 when only optimizing for alignment. 

It is interesting to note that the maximum alignment obtained is the same as in the first case, with (cos^ 0) (tf) = 
0.867 fFig. llOfl . even though the two wave functions obtained involve very different mixtures of spherical harmonics. 
One striking contrast between the two is actually not visible when looking only at (cos^ O): in the second case, the 
state obtained is strictly aligned, in the sense that the angular distribution is symmetric with respect to 6 — 7r/2. 
In the first case, the maximum in alignment corresponds also to a maximum in orientation, with {cos9) (if) = 0.891, 
whereas in the second case (cos0) (if) = —0.027. 



C. Considering a rotational temperature 



By starting all simulations from the ground rotational state, we have in fact made the approximation of a zero initial 
rotational temperature. From previous work on laser-induced alignment and orientationfi^s^iSSi^iiS it is known that 
considering a higher, experimentally more realistic initial rotational temperature will lead to an important decrease in 
the amount of orientation/alignment obtained. This can be seen in Fig. 1111 where we show how orientation is affected 
when the optimal field presented in Fig. ^a) is applied to a thermal ensemble of initial temperature k-oT/B « 4.77 
(corresponding to 10 K for HCN). 

This difficulty can be overcome by adapting the optimization method to take into account the initial thermal 
distributioni^SiS We present here the optimization of orientation, which can be trivially extended to the case of 
alignment or of a suitable mixed-state target, 52. The measure of orientation now reads 

oo 

((cos0))(t) = Q-i^exp 
j=o 



-BJ(J+\) 



^ (cos6i)j^A,^ (i). 



(14) 



where 



Q = ^(2J+ l)exp 



,7=0 



-BJ(J+ 1) 



knT 



(15) 



is the partition function and {cos9)j (t) = {^pJ,M{0'^t)\ cosd IV'j.m (^; i)) is obtained by solving the TDSE with the 
Hamiltonian for the initial condition V'j,m(^; t = 0) = Yj,m((^)- Considering a temperature k-QT/B « 4.77, we can 
restrict the sum over J in Eq. H14|) to Jmax = 7, with 16 basis functions [see Eq. 0] used for the time evolution. The 
monotonic algorithm (|10|) becomes 



£^+^{t) = -Q-i£cxp 



,7=0 



-B.7(.7+l) 



^5+J(0,i = O) 



B - ^iQ£''+^{t) cos6'- ^ 

Y.,,M{e). 



J 

V - 

M=-J 



?(V.;+J(t)|A.o cose|x",,M(*)) 



Aa cos^ 6 + 



'yj.M 



x5,M(t))-A(t) 



(16) 



The resulting optimized field, given in Fig. I12f a). allows to reach {{cos9)) (if) — 0.686, which is much better than 
the value of 0.434 obtained when the T — K field is applied to the same thermal distribution, as shown in Fig. ^2 
Looking at the Fourier transform of the field. Fig. I12r b). we see that the mechanism to achieve orientation is slightly 
different, with the lower frequencies present for a much longer time than previously observed, which reflects the fact 
that many J ^ states are initially populated. Nevertheless, the ladder climbing structure of the resonances is still 
present. 



IV. CONCLUSION 



Using a monotonically convergent algorithm, we have searched for the optimal electric field maximizing cither the 
alignment or the orientation of a linear molecule, taken in a rigid rotor approximation. We have carried out the 
optimization both in terms of maximization of observables corresponding to orientation/alignment and using target 
states offering a good compromise between the efficiency of orientation/alignment and its duration. 

We have found that, starting from the ground rotational state, the optimal fields allow to reach orienta- 
tion/alignment by rotational ladder climbing, i.e., by successive resonant excitation of neighboring rotational levels. 
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We insist on the fact that this scenario appears "naturally" from the physics of the problem and is not imposed a 
priori by the optimization algorithm. This process allows to reach an orientation of (cos0) = 0.909 or an alignment 
of (cos'^ 6*) = 0.867. Target states can also be reached to within better than 0.5%. 

We have also shown how our optimization method can be extended to the more realistic case of an initial thermal 
distribution of rotational states. Using orientation as an illustrative example, we obtained a value of {{cosO)) = 0.686 
at a temperature k^T / B « 4.77 (corresponding to 10 K for HCN). This result is much better than those previously 
reported,-^ even though we are considering here a higher rotational temperature. 

As a rigid rotor model was used, this study did not take into account any vibrational excitation, which could hinder 
or enhance the orientation/alignment obtained. By including vibrations into the model, it is possible to use different 
control paths not involving direct rotational excitation, enabling the choice of infra-red lasers as control fieldsiSiiSiS 
Vibration-rotation coupling can also lead to cross-revivals of vibrational wave packetsiSSiSi Future work will thus take 
into account the vibration of the molecule. 
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FIG. 1: (a) Electric field obtained with criterion jTi for the optimization of (cosO). (b) Short-time Fourier transform of the 
field in (a). 



9 




FIG. 3: Orientation, as measured by {cos 9), obtained for a rigid rotor interacting with the electric field given in Fig.^a). The 
field-free evolution is then periodic with period Trot- 




FIG. 5: Orientation, as measured by {cos 6), obtained for the interaction with the electric field given in Fig. ^a) (solid line) 
and Fig. 0Ja) (dashed line). The field- free evolution is then periodic with period Trot. 




FIG. 7: Time evolution of the population of rotational states of a rigid rotor interacting with the electric field given 
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FIG. 11: Orientation, as measured by ((cosS)), obtained for a rigid rotor interacting with the electric field given in Fig. 
starting from the rotational ground state (solid line) and from a rotational temperature of ksT/B ~ 4.77 (dotted line). For 
comparison, the orientation obtained with the field optimized for this temperature [Fig. I12f a'l1 is also shown (dashed line). 




FIG. 12: (a) Electric field obtained for the optimization of ((cos9)) starting from a rotational temperature of ksT/B « 4.77. 
(b) Short-time Fourier transform of the field in (a). 



